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Abstract. We apply the pixel method of analysis (some- 
times called "pixel lensing") to a small subset of the 
EROS-1 microlensing observations of the bar of the Large 
Magellanic Cloud (LMC). The pixel method is designed 
, to find microlensing events of unresolved source stars and 
had heretofore been applied only to M31 where essentially 
, all sources are unresolved. With our analysis optimised 
for the detection of long-duration microlensing events due 
to 0.01-1 Mq Machos, we detect no microlensing events 
and compute the corresponding detection efficiencies. We 
show that the pixel method, applied to crowded fields, 
should detect 10 to 20 times more microlensing events for 
M > O.O5M0 Machos compared to a classical analysis of 
the same data which latter monitors only resolved stars. 
In particular, we show that for a full halo of Machos in 
the mass range 0.1 - 0.5 Mq, a pixel analysis of the three- 
year EROS-1 data set covering 0.39 deg 2 would yield ~ 4 
events. 

Key words: Methods: data analysis - Techniques: pho- 
tometric - Galaxy: halo - Galaxies: Magellanic Clouds - 
Cosmology: dark matter - Cosmology: gravitational lens- 
ing 
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1. Introduction 

Microlensing searches can probe the distribution of MAs- 
sive Compact Halo Objects (Machos) in the dark halo 
of our Galaxy or more distant galaxies (Paczyhski 
(1986,1996), Griest (1991)). Several lines of sight are now 
under investigation, and events have been claimed in sev- 
eral directions: towards the LMC (Alcock et al. 1997a, 
Renault et al. 1997), the SMC (Alcock et al. 1997b, Palan- 
que-Delabrouille et al. 1998), in the direction of the Galac- 
tic Bulge (Alard et al. 1997, Alcock et al. 1997c, Udalski 
et al 1994) and more recently towards spiral arms (Derue 
et al. 1998), giving some first evidences of the Macho dis- 
tribution towards these lines of sight. These results are 
based on a star monitoring analysis: the fluxes of several 
millions of resolved stars are monitored. As first discussed 
by Crotts (1992) and Baillon et al. (1993), events due to 
unresolved stars essentially escape detection of these anal- 
yses. Such stars, beyond the crowding limit or too dim 
to resolve, could significantly contribute to the number of 
detectable events. This is illustrated by the detection of 
two variable objects in the MACHO analysis (Alcock et al. 
1997a) of LMC data, which could not be resolved at their 
minimum luminosity. The detection of the variation was 
nevertheless possible because the reference images used 
to establish the catalogue of monitored stars were taken 
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during their maximum luminosity, when the stars were re- 
solved. However, such events can neither be kept for fur- 
ther considerations nor be included in the computation of 
detection efficiencies, in the star monitoring approach. 

In this paper, we apply a pixel analysis to the 
EROS 91-92 data (10% of the whole EROS-1 CCD data 
set) of the LMC Bar. The idea is to monitor the flux of 
all the pixels present on the images, thus achieving a good 
sensitivity to the whole stellar content of the images. The 
magnification of one unresolved star can be detected as 
a variation of the pixel flux, provided that the magnifi- 
cation is high enough. - In the following, we will refer 
to this approach as pixel monitoring, as opposed to star 
monitoring referring to classical analyses restricted to re- 
solved stars. - The main uncertainties of this approach 
concern our ability to account properly for variations of 
the observational conditions, and to be able to disentan- 
gle intrinsic variations from observational systematics. In 
Paper I (Melchior et al. 1998a), dedicated to the descrip- 
tion of the treatment of the data and the production of 
2.1 x 10 6 super-pixel light curves, we have shown that an 
average stability of 1.8% of the super-pixel flux is achieved 
in blue and 1.3% in red, about twice the expected photon 
noise. This homogeneous set of super-pixel light curves 
is called AGAPEROS: each of these light curves covers 
a period of 120 days and is composed of about 90 mea- 
surements. With this rather short period of observation, 
we show here how it is possible to investigate the Macho 
mass range of interest (M ~ 10 -2 — 1M Q ) with the ex- 
isting EROS-1 CCD data set, initially designed for short 
time scale events in a mass range (M ~ 10~ 8 — 1O _3 M0) 
where no event has been detected. 

Microlensing selection with the Pixel Method 

We first present the simple formalism used to describe the 
pixel events in which we are interested. The pixel flux cj) p , 
affected by a microlensing event, can be written as: 



the Einstein radius: 



<% = OiAiuW + (p. 



(1) 



where i is the measurement number, t is the time, a is 
the seeing fraction, A the magnification, <fi* the flux of 
the star of interest at rest (i.e., unmagnified) , and (f> hs 
includes the sky and stellar backgrounds. The magnifica- 
tion (Paczyhski 1986) depends on the normalized impact 
parameter u(U): 



with 
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where vt is the Macho transverse velocity, Uq and to the 
impact parameter and time of closest approach, and Re 
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The typical time scale of the variation is the Einstein ra- 
dius crossing time: 

t E = ^ (3) 
vt 

In Paper I, we showed that the variations of the ob- 
servational conditions, which obviously affect Eq. [jj can 
be corrected: each image is first geometrically aligned with 
the reference image. Then sky background and absorption 
factor are corrected to the values of the reference image. 
Finally, each super-pixel flux is adjusted to take account of 
the seeing variation, affecting on. Since the mean seeing is 
about 3 arc-second, we estimate a to be on average 0.7 for 
the corrected 3"6 x 3"6 super-pixel light curves, obtained 
in Paper I. We showed that in the absence of microlensing 
events (^4 (ti) = 1), we achieve a proper understanding of 
the errors affecting these light curves, and that Eq. || de- 
scribes to a good approximation the light curves we are 
studying. 

The usual requirements used for the selection of mi- 
crolensing events detected by star monitoring can be ap- 
plied here: 

— As the microlensing phenomenon is a transient and 
rare phenomenon, it should produce a unique signifi- 
cant variation in the star flux. 

— It must be achromatic. This characteristic has two ap- 
plications for a pixel analysis: the time of the maximum 
has to be the same in both colours and the ratio 



MB) 



bg 



(4) 



must remain constant, during the variation. 

— The shape must be compatible with Eq. |l|. 

The first criterion allows us to remove recurrent variable 
stars as well as most of the noise, while the two other 
criteria will be applied to the few remaining light curves 
at the end of the selection process. 

— Last, we consider the fact that the probability for a 
star to be lensed is independent of its type. This will 
allow us to reject specific populations of variable stars. 

In Sect. |^, we present the Monte-Carlo simulations 
used in this article. Then, in Sect. |L we describe step by 
step the selection procedure designed to detect microlens- 
ing events and applied to these light curves. In Sect. § 
we show how the selected variations can be eliminated as 
compatible with variable stars. In Sect. ^ we discuss the 
detection efficiencies achieved by this analysis, and com- 
pare the number of expected events with the sensitivity 
of star monitoring. We rely on these results in Sect. ^| to 
discuss the possible prospects of this approach. 
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2. A useful tool: mock super-pixel light curves 
with microlensing events 

Monte-Carlo simulations described in Baillon et al. (1993) 
gave a first estimate of the number of events expected with 
a pixel analysis. The main uncertainties discussed there 



derived from the noise present in real data. In Sect. 2.1 



we present a summary of the simulations of microlensing 
events used in this paper. These model the physical ingre- 
dients including the halo density profile, the luminosity 
function of source stars, and the Macho velocity distribu- 
tion and mass function (see Baillon et al. 1993). In Sect. 
2.2, we define a minimal threshold that will be useful later 
on for the interpretation of the results of our analysis. In 
Sect. 2.3, we discuss the characteristics of the simulated 



events as expected for an ideal analysis of our data set. 
Last, in Sect. 2.4, we add to this model the characteris- 



tics of the AGAPEROS data, and thus produce realistic 
mock super-pixel light curves. This tool will be used in 
the following to adjust the selection criteria in Sect. ||and 
to compute the detection efficiencies in Sect. [| 



2.1. Physical ingredients 

We assume an isothermal halo with a core radius of 5 
kpc, normalized at the solar neighbourhood to pq = 
8. x l(T 3 Af©pc- 3 (Flores 1988) and filled with compact 
objects with a given mass M as discussed by Griest 
(1991). We adopt a LMC distance of 50 kpc. The cor- 
responding optical depth for a full halo is the same as 
the one used by the MACHO group (Alcock et al. 1997a) 
Tfuii = 4.7 x 10~ 7 . It is to be noted that our estimate of 
the expected number of events assumes a full halo and it 
should be multiplied by a factor / < 0.5 according to the 
MACHO and EROS results (Alcock et al. 1997a, Renault 
et al. 1997), where / is the halo fraction actually filled 
with Machos. According to the MACHO results (Alcock 
et al. 1997a), it is most probably smaller than 0.5. Note 
that no halo flattening has been considered at this stage, 
but more sophisticated models could be implemented and 
tested once serious candidates are detected. 

We calculate the number of potential lenses with a 
fixed mass M located in the cone pointing towards our 
field of view. We assign a random position to the Ma- 
cho and choose its transverse velocity vt from a two- 
dimensioned Maxwellian distribution g(vx) weighted by 
vt- For each given Macho, we determine the probability 
that a star will lie close enough to this line of sight to 
give rise to a microlensing event. We use the luminosity 
function described in Baillon et al. (1993) which is based 
on Hardy et al. (1984) for the bright stars, on Ardeberg 
et al. (1985) up to magnitude V — 23 and is extrapo- 
lated to the faint end using the luminosity function of 
the solar neighbourhood (Allen 1973). Note that on the 
one hand, the details of this latter extrapolation are not 
important because, as we show below, few events are de- 




Fig. 1. Luminosity function used for the simulations: the 
number of stars per arcsec 2 or per pixel ( 1 ." 2 1 x l."21), 
normalised to a surface brightness /iy = 21 is exhibited 
as a function of the V magnitude. This shows the typical 
stellar content of a pixel. 



tectable for sources fainter than V — 24; but on the other 
hand, the connection between these 3 sets of observations 
is a source of uncertainties. This luminosity function is dis- 
played in Fig. |l| and is quite compatible with recent mea- 
surements (Holtzman et al. 1997, Ardeberg et al. 1997). 
We then normalize this function to a surface brightness 
of /iy = 21 (de Vaucouleurs 1957). The star's magnitude 
is drawn from a uniform distribution and the simulated 
event is then weighted according to the luminosity func- 
tion. As described in Baillon et al. (1993), we account for 
possible finite source effects that are expected when the 
stellar radius projected onto the plane of the Macho is 
comparable with the Einstein radius. We are then able to 
compute the number of expected events using well-known 
Monte-Carlo integration techniques. 

2.2. Minimal requirements for simulated microlensing 
events 

We require here some minimal requirements that will de- 
fine a set of simulated microlensing events that could be 
detected with an ideal experiment. 

As we do not expect to detect a significant num- 
ber of low-magnification events, we introduce a thresh- 
old A max > 1.34 in our simulations. Detectable low- 
magnification events affect bright stars and hence would 
have already been detected by the previous EROS star 
monitoring analysis anyway. Moreover, the magnification 
by such a small factor of a dim star would be completely 
buried into the noise, therefore one must add a visibility 
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Fig. 2. Initial set of the simulated light curves with some 
microlensing events. The distributions of the impact pa- 
rameters uq (panel a) and of the V magnitude of the un- 
magnified star (panel b) are exhibited. The histograms 
shown are computed for 0.5 Mq Machos. At this stage 
simulations with different Macho's mass will give the same 
histograms but with a different normalization. 



condition. We choose the following: at the time of max- 
imum magnification, the flux of the central super-pixel 
of a magnified star should rise higher than 3<r above the 
background, a being taken as twice the photon noise. It 
is important to note that this threshold does not depend 
on the duration of the event. It only removes events that 
we would not detect in any case. 

The effect of this requirement on the characteristics of 
the simulated sample can be seen in Fig. ||. The impact 
parameter distribution is no longer expected flat. This 3<r 
threshold introduces a (necessary) bias into the impact 
parameter distribution towards small values. The majority 
of the events are expected to affect dim stars with a small 
impact parameter u < 0.2. 

The simulations including these two thresholds (A > 
1.34, and S/N > 3cr) will be used as a reference for the 
computation of our detection efficiency in Sect. [|. 

2.3. Preliminary estimates 

Firstly, we discuss the number of microlensing that can be 
expected with an ideal experiment. Secondly, we estimate 
the number of stars effectively monitored in our reference 
set. 



Table 1. Number of expected microlensing events as a 
function of the Macho mass, given the minimal require- 
ments discussed in the text. 



M/Mq 


0.01 


0.05 


0.1 


0.5 


1.0 


Ar AGAPEROS // 


9.8 


4.4 


3.1 


1.4 


1.0 



2.3.1. Number of expected events 

The number of microlensing events estimated by our sim- 
ulations for a halo filled with Machos of mass M is thus 
defined as: 



JVcvt = / d<j) 



dD 



Altars (</>*) X T ohs / dv 



vrgivr) x 



2u (<j>*)R E (D) 



f P(D) 
M 



,t E (M,D,v T )) 



(5) 



where N silirs ((j)*) is the true number of stars with a flux 
between (jf and 4>* + d<jf present in the sky area studied, 
T bs is the duration of the observations (120 days), p(D) 
the Macho density distribution, D the position of the Ma- 
cho and e is the efficiency (e = 1 for an ideal experiment) . 

Tab. gives the number of microlensing events that 
can be expected with an ideal analysis of the super-pixel 
light curves produced in Paper I. In the large-mass range, 
where all the known microlensing candidates on the LMC 
have been identified, we expect between 1 and 10 mi- 
crolensing events assuming a full halo. 

2.3.2. Effective number of monitored stars 

If we were able to select all light curves of our reference 
sample, we would define the equivalent number of moni- 
tored stars A^agaperos as follows: 



N. 



AGAPEROS _ 
stars 



#^Btara(^*)«0(^*|3(T) 



(6) 



U(j((j)*\3a) is the threshold impact parameter that enters 
Eq. [| which accounts for the 3a deviation imposed at the 
time of maximum magnification in Sect. [2.2] . Actually, this 
number A^^ PEROS would only depend on the luminos- 
ity function and the definition of our reference sample. 
Hence, we can consider that we effectively monitor the 
equivalent of 2.2 x 10 6 stars, whose mean magnitude is 
22.1 (see Fig. |b). 

If one integrates the luminosity function]^] (Fig. over 
a pixel area ( 1 2 1 x l."21), one finds one star in the mag- 
nitude range 20 — 24, that could undergo a microlensing 



We estimate that the remaining uncertainty on the LF 
translates into 20% uncertainty in this effective number of 
stars. 
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Fig. 3. Two examples of simulated super-pixel light curves 
containing a microlensing event. 

event. This explains why the effective number of stars thus 
defined is of the same order as the number of pixels. 

24. Model of the data 

The idea is to simulate super-pixel light curves that in- 
clude microlensing events. We compute the flux of the star 
- affected by a microlensing variation - which enters the 
3."6 x 3."6 super-pixel. We also add a background flux - 
together with expected read-out and photon noises - in or- 
der to obtain realistic mock light curves. The computation 
of these fluxes takes account of the pass-band of the filters, 
the quantum efficiency of the CCD camera and its gain. 
(See Paper I and references therein for more quantitative 
information about the characteristics of the raw data.) Ac- 
tual spacing and variations of the observational conditions 
(absorption and sky background), measured on the data, 
are also taken into account in this procedure. Similarly 
to what is done for real data, the averaging procedure of 
the measurements available each night is applied to these 
mock light curves, as well as to the computation of error 
bars. We multiply these errors, assumed to be gaussian 
and uncorrelated, by the factor found in Paper I between 
the measured dispersion and the expected photon noise. 

We finally get mock light curves typical of the mi- 
crolensing events we are looking for. Fig. [5] displays two 
examples of typical expected events. 



3. Selection of microlensing events 

We apply a pixel analysis designed to select microlensing 
events on the EROS 91-92 data of the LMC. Using the 
methods described in Paper I and applied on the AGAP- 
EROS data set, we constructed some 2.1 x 10 6 (real) super- 
pixel light curves, cleaned of all observational effects and 
corrected for systematic effects to the degree possible. In 
this section, we define a selection process designed to se- 
lect possible microlensing events. The application of a ba- 
sic trigger - detection of bumps - reveals a large number of 
variations, most of them corresponding to obvious variable 
stars, but also to some noisy variations that we will have 
to eliminate. Owing to the averaging of the images taken 
each night (see Paper I), some of the variations are most 
probably due to short-time scale variables already studied 
elsewhere (Beaulieu 1995, Grison 1995). As Renault et 
al. (1997,1998) have already excluded the small-mass Ma- 
cho range with this data set, we choose to optimise our 
sensitivity to long time-scale variations, corresponding to 
the large-mass range in which all the known microlensing 
events have been detected. 

The selection criteria should remove intrinsic varia- 
tions and systematic effects, keeping genuine microlens- 
ing events. We describe these various criteria successively 
applied to the data. Some rather loose cuts, applied on 
super-pixel light curves, turn out to be sufficient to re- 
duce considerably the number of light curves to analyse. 
Then a visual inspection of the 120 remaining curves con- 
firms that they are affected by genuine variations and they 
will be further studied in Sect. The efficiency of our se- 
lection procedure with simulated super-pixel light curves 
has been checked in both colours. 

We work on a simulation based on 5 x 10 6 realizations, 
which allows to have small statistical errors on the num- 
ber of expected events. The number of events given at 
each step of the selection procedure is calculated assum- 
ing a halo full of 0.5M Q Machos. In Sect. [| we discuss 
the sensitivity actually achieved in a larger-mass range 
0.01M Q < M < 1M Q . 

The selection procedure, defined in this section and 
summarised in Tab. |^, splits into 3 steps: (1) a signifi- 
cant variation must be present in at least one colour. (2) 
Light curves with any significant secondary variations are 
eliminated. (3) A correlation between the two colours is 
required. The two first thresholds are based on the follow- 
ing definition of the variations or "bumps" . 



Definition of a bump A baseline <pu is calculated for each 
super-pixel light curve as the minimum of a running av- 
erage over 5 successive flux measurements. o~bi is the error 
associated with the baseline flux determination. All the 
light curves are scanned for the detection of bumps, de- 
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Table 2. Summary of the selection procedure: the per- 
centage of events kept at each steps is given for the simu- 
lations (2nd column) and the data (4rd column). The 3rd 
and 5th columns give the number of simulated and real 
light curves kept. 



EROS - CCD no03 - pixel (045,400) 



Starting from 



Simulations 

1.4 

events 



Data 

2.1 x 10 6 
light curves 



One bump 

-Li> 500 51.0% 0.71 0.25% 5 338 

- 3 points above 3a 71.2% 0.51 52.3% 2 789 

No second bump 90.3% 0.46 77.5% 2162 

Good correlation 85.5% 0.39 5.5% 120 



fined as at least 3 consecutive measurements lying above 
the baseline by at least 3a n 



<Jbr 



(7) 



where a' n is the error associated to each super-pixel flux 
computed in paper I for the night n. The bump ends when 
at least 2 consecutive measurements lie below this thresh- 
old. Each bump i is characterised using a likelihood func- 
tion: 



U = - In 



( II P&>4>n) given j 



(8) 



4> n is the super-pixel flux for the measurement n, whose 
computation for all the light curves is detailed in Paper I. 

3.1. At least one significant variation 

Firstly, we require one significant bump in at least one 
colour. Secondly, we then look for a minimal variation in 
the other colour. 



A large bump in at least one colour In order to identify 
significant variations, we require the likelihood function L\ 
associated with the largest fluctuation to be larger than 
500 for at least one colour. This value is chosen using the 
Monte-Carlo simulations to optimise the S/N ratio. In a 
given colour, we search for clusters of super-pixels^] hav- 
ing each L\ larger than 500. We then select the central 
super-pixel of each cluster, if it also satisfies the L\ > 500 
requirement. The latter requirement is intended to remove 
some artifacts, in particular close to bright stars. The next 
cuts will be applied on these central super-pixels. 



We use a Friend of Friends algorithm (see for instance 
Huchra & Geller (1982)). 
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Fig. 4. Super-pixel light curve with several variations re- 
jected with the cut on L 2 (if = 1366. ; Lf = 2958. ; 

= 565. ; L 2 = 1655.). The upper curve exhibits the 
light curve in blue and the lower one in red. 



Minimal variation in the other colour The previous cut 
has allowed to detect significant variations in at least one 
colour. We now require that the first fluctuation of the 
central super-pixels of the other colour to have at least 3 
consecutive points 3tr above the baseline. Although it cor- 
responds to a quite small value of L (~ 15), it constitutes 
a first requirement of achromaticity. 

3.2. No significant second bump 

At this stage, we have selected super-pixel light curves 
with at least one significant variation. Now, it is impor- 
tant to check uniqueness. We then require the second most 
significant fluctuation to have L2 < 250 in both colour. 
Fig. ^ shows an example of rejected light curves with two 
variations. 

3.3. Correlation between the two colours 

In order to select the expected long time-scale variations, 
we choose for the next requirement a reasonable correla- 
tion between the blue and red light curves. 

Given the constraints mentioned above, the correlation 
Pcoi between the two colours: 



v E„(^-(0» 2 bluc E n (^-(^>f 



(9) 



achieves a good sensitivity to the achromaticity and to the 
dispersion of the measurements. As shown by the mock 
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Fig. 5. Histogram of the correlation factor computed be- 
tween the two colours. The hatched area corresponds to 
super-pixel light curves with a correlation factor below 
the threshold p = 0.8. The histogram superimposed with 
a dashed line on this figure corresponds to the simula- 
tions and scales with the tick marks and values given on 
the right axis. 



light curves in Fig. || (dashed curve), a clear correlation 
between the two colours is expected: most of the mock 
curves (85.5%) lie above a threshold p = 0.8. Fig. || dis- 
plays the corresponding histogram for the real light curves 
selected so far (full line). This distribution is quite differ- 
ent and exhibits a peak close to p = 0. With a threshold 
at p^ ol > 0.8, 94.4% of the remaining light curves are re- 
moved: 120 variations remain. Fig. ^ displays one of these 
light curves: its shape is in quite good agreement with 
what we can expect from a microlensing event although 
the period of observation is short compared to the dura- 
tion of the variation, and it is not possible to test up to 
now the uniqueness of the variation. Fig. ^ displays an- 
other light curve, one whose shape is clearly incompatible 
with a standard microlensing event. Fig. || shows a light 
curve for which it is impossible to draw any conclusion 
based only on compatibility with the microlensing shape: 
the period of observation is much shorter than the time- 
scale of the variation. 

When the criteria described above are applied in the 
simulations, we expect 0.38 x / events for 0.5M Q Machos 
filling a fraction / of the halo. It is obvious then that the 
120 selected light curves are in clear excess with respect 
to what is expected and need to be further studied. 



EROS - CCD no03 - pixel (289,296) 

blue 



1200 



1 100 



ft 



iiiiU*^ 1 



20 40 60 80 100 120 

JD-24486608.5 



-e- 8000 f 



red 



7000 - 



T1V 



6000 



20 40 60 80 100 120 

JD-24486608.5 

Fig. 6. Super-pixel light curve with a shape compatible 

with a microlensing event. (if = 641. ; if = 744. ; Lf = 
R 

2 



57. ; Lf = 0. ; p = 



EROS - CCD no03 - pixel (059,180) 



: blue 

1000 - 



950 



IT 



850 Hi!" n\ 

20 40 60 80 100 120 

JD-24486608.5 



-e- 






7000 


: red 




6000 




i y f ) 








5000 




J" / 


4000 






3000 








20 40 60 80 100 120 

JD-24486608.5 

Fig. 7. Super-pixel light curve with a shape incompatible 
with a microlensing event (if = 1386. ; Lf = 16416. ; 
if = 80. ; if = 0. p = 0.93) 



4. Colour magnitude diagram 

Since the expected background for microlensing events is 
due to variable stars, we have first to determine the posi- 
tion of the 120 light curves with variations on the colour 
magnitude diagram (CMD), before drawing any conclu- 
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Fig. 8. Super-pixel light curve for which it is inconclusive. 
(Lf = 1791. ; if = 12927. ; Lf = 0. ; Lf = 176. p = 
0.96) 

sion on the nature of these variations. The problem is now 
to know how we can estimate the magnitude of the under- 
lying star responsible for the observed variation. 

4--1- Magnitude determination 

Using an image with a seeing close to the mean value 
(3 arc-second) - for which no significant seeing correc- 
tion is required - we perform photometry on the sur- 
rounding stars using DAOPHOT which returns the total 
background lying within 2 arc-second of the center of the 
super-pixel. We then subtract this background from the 
super-pixel flux measured when the star is at maximum 
luminosity, and account for the seeing fraction of the star 
flux entering the super-pixel. Hence, we deduce the mag- 
nitude of the star at the maximum. Due to the crowding 
of the LMC Bar, this aperture photometry is the most effi- 
cient way to estimate a magnitude for the stars responsible 
of the detected variations. More details about magnitude 
determination will be further addressed in Melchior et al. 
(1998b). 

This estimate is mainly intended to study the position 
in the CMD of the dominant source of flux of the varying 
pixel and in particular if it lies in marginal locations of 
the colour-magnitude diagram, characteristic of variable 
stars. We estimate the uncertainties on this magnitude 
determination as the square root of the sum of the squares 
of the two following components. The first is the error on 
the super-pixel flux. The second one is estimated as 10% 
of the "star" flux and is expected due to uncertainties in 
the star position inside the inner pixel of the super-pixel. 
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Fig. 9. Colour magnitude diagram: the small dots corre- 
sponds to the stars detected by the EROS-1 star monitor- 
ing analysis (Renault, 1996). The circles corresponds to 
the 120 selected variations discussed in the text. Among 
them the filled circles are Miras detected by Hughes (1989) 
in the I band (see text). 

Table 3. Miras: matches with a previous study by Hughes 
(1989). The following characteristics are given: the posi- 
tion at 1950.0 epoch, the CCD number, the pixel position, 
the mean / magnitude, the period in days 
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In extreme cases - when the star flux at maximum is not 
the main contribution of the super-pixel flux - these errors 
can be underestimated. 



4-.2. Discussion 

Fig. ^| displays a CMD with the stars detected by the 
EROS group (dots) and the underlying stars (circles) as- 
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Fig. 11. Map of the surroundings of 
the Mira whose light curve is dis- 
played on Fig. [h]: each figure dis- 
plays the flux as a function of posi- 
tion on an area of 20.6" x 20.6". Up- 
per figures a and b correspond to 
blue images, and lower ones c and 
d to red images. Figures a and c 
on the left side exhibit the star at 
the maximum of the variations. On 
the figures on the right side b and 
d, this star has disappeared in the 
background. Note that such an event 
would have escaped any star moni- 
toring, and could not have been ac- 
counted for as part of a blend. 



sociated with the selected super-pixel light curves. The 
variations kept by our analysis are not representative of 
the bulk of the stars: all of them but two lie in the red 
part of the CMD. 



The 118 red variable stars The red detected variable 
stars are all located in the same area of the colour mag- 
nitude diagram. We have even been able to check that 
10 of them have already been recorded by Hughes (1989) 
as Miras with a study in the I band. They are displayed 
with filled circles on Fig |] and their known characteristics 
are displayed in Tab. [| Fig. |l^ exhibits the light curve of 
one of these Miras, with a 163 days period. Not surpris- 
ingly the variation is significant but only sampled over 120 
days, and this is the shortest period of the Miras listed in 
Tab. |[ Two of these Miras have been shown previously. 
One is presented in Fig. 0, it has a period of 650 days. An- 
other one is exhibited in Fig || with a period of 453 days. 
The catalogue of Hughes (1989) contains 41 variable stars 
overlapping the field studied here and most of them have 
been rejected at an earlier stage of this analysis. The thick 
lines show the area of the CMD excluded by the EROS 
group corresponding to the regions where variable stars 
are expected. It is then highly probable that the other red 
variable stars selected are also Long Period Variables, as 



they lie in the same area of the CMD. The fact that the 
majority of the red variable stars selected by this analy- 
sis have not been identified previously demonstrates the 
potential interest of the pixel method for the detection of 
Long Period Variable (LPV) stars with respect to clas- 
sical analysis restricted to the study of resolved stars. A 
comprehensive analysis of these variable stars rejected as 
background of the microlensing search will be presented 
elsewhere (Melchior et al. 1998b). 

About 10 of these red variable stars lie below the 
crowding limit and Fig. [lO] displays the light curve of one 
of them: the Mira already discussed above. For these stars, 
we are not able to detect variations around their minimum 
flux. This is illustrated in Fig. [ll] which shows the field 
surrounding this star: an unresolved star has exhibited a 
variation. Although this particular example would have 
escaped the EROS-1 star monitoring applied on the same 
data, it has already been identified in another wavelength 
(I band) by Hughes (1989). 

The 2 blue variable stars The brightest blue variable 
star (with Be — 15) shown in Fig. [l^ belongs to the 
sample of pre-main-sequence stars selected by Beaulieu 
et al. (1996). The short time scale variation which super- 
imposes on long-time scale 0.3 magnitude variation is real 
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Fig. 10. Super-pixel light curve with a variation detected 
by our analysis which has been identify as a Mira. (Lf = 
10188. ; Lf = 47940. ; Lf = 112. ; Lf = 0. p = 0.96) 

and this feature excludes the simple microlensing inter- 
pretation anyway. 

The other blue variable star is characterised by a small 
amplitude (ABe < 0.3) and a duration longer than 120 
days. These features together with the position of this 
variable in the CMD indicate that the variation is com- 
patible with the new class of variable stars called the blue 
bumpers identified by the MACHO group (Cook et al. 
1995, Alcock et al. 1996). 

4-3. Final cut on the colour magnitude diagram 

These considerations provide convincing evidences that 
the selected variations are variable stars (Long Period 
Variables for most of them). We decide to apply the same 
cut as EROS on the CMD, displayed with the thick lines 
shown in Fig. ||. These regions are filled with a negligible 
number of stars (< 1.3%), that are moreover expected to 
be variable. The elimination of this area does not signifi- 
cantly affect our sensitivity to microlensing events which 
are expected to occur independently of the star's position 
in the colour-magnitude diagram. 

5. Results 

Although we detect no microlensing events with the selec- 
tion procedure described above, we would have detected 
them if there were some. Whereas the detection of variable 
stars gives a first idea of the sensitivity achieved by this 
analysis, the Monte-Carlo simulations provide detection 
efficiencies. Firstly, we present the detection efficiencies 
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Fig. 12. Super-pixel light curve with a variation affecting 
a blue star detected by our analysis. It has been detected 
previously by Beaulieu et al. (1996) and identified as a 
pre-main-sequcnce star candidate labelled ELHC2. (Lf — 
1747. ; Lf = 3850. ; if = 0. ; = 146. p = 0.89) 

Table 4. Percentage of events kept for each step of the 
selection procedure for different Macho's mass. The last 
row provides the percentage of events kept with respect 
to the initial set of events. 
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Li > 500 in at 


32.0% 


42.1% 


43.7% 
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48.4% 


least one colour 












3 pts above 3a 


44.1% 


58.6% 


65.0% 


71.2% 
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in both colours 












L 2 < 250 in 


96.4% 


93.3% 


92.9% 


89.7% 


89.4% 


both colours 












p > 0.8 


75.9% 


76.0% 


75.0% 


85.5% 


88.7% 



Total efficiency 10.3% 17.4% 19.8% 26.7% 27.9% 



achieved for this pixel analysis. Secondly, we compare our 
sensitivity with those achieved by star monitoring analy- 
ses. 

5.1. Detection efficiencies for pixel monitoring 

In Sect. |^, we detail the effect of the selection procedure 
on simulated light curves. For the clarity of the discus- 
sion, we restricted then the comparison to events expected 
with a halo full of 0.5M© Macho. We study here the sen- 
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Fig. 13. Detection efficiencies computed as a function of 
the Einstein radius crossing time tE- 




26 28 



Fig. 14. Selected simulated light curves. The distributions 
of the impact parameter uo (a) and of the V magnitude 
(b) of the star at rest are exhibited, for 0.5 Mq Machos. 
Normalisation corresponds to a full halo. 




Fig. 15. Selected simulated light curves. The distributions 
of the impact parameter uq (a) and of the V magnitude 
(b) of the star at rest are exhibited, for 0.01 Mq Machos. 
Normalisation corresponds to a full halo. 
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Fig. 16. Detection efficiencies computed as a function of 
the magnitude of the un-magnified star. 



sitivity achieved by our selection procedure on a larger- 
mass range. Tab. || shows the percentage of simulated 
events selected by each step of our selection procedure for 
0.01M© < M < IMq. It appears clearly that we have op- 
timised the selection procedure for this Macho mass range, 
and that efficiency is lost with decreasing Macho mass. 
The last cut on the correlation factor (p > 0.8) is less 
favourable for M < 0.1M Q : for variations affecting only 
part of the period of observation, it is more optimal to 
restrict the computation of this coefficient to the portions 
of the light curve undergoing variations. These efficiencies 
can be studied as a function of the duration tE as shown 
in Fig. |l^. Due to the temporal sampling, the selection 
procedure is less efficient for short duration events. The 
efficiency remains constant for long-duration events, as we 
do not require a stable baseline. 



Fig. [lj shows the distribution of impact parameters 
and V magnitude for simulated events due to 0.5 Mq Ma- 
chos satisfying all our requirements. It is clear that with 
respect to Fig. ^J, the selection procedure eliminates mi- 
crolensing events affecting very dim stars (V > 24), and 
that the main contribution is expected due to events af- 
fecting dim stars with a small impact parameter. Fig. ^ 
gives the same histograms but for events due to 0.01M Q 
Machos. Events affecting dim stars are much more diffi- 
cult to detect if they are short. Hence, detected events 
affect on average brighter stars and the impact parameter 
distribution appears flatter. 

Fig. |l6| exhibits the detection efficiencies achieved as a 
function of the V magnitude of the un-magnified star and 
illustrates the previous point. The longer the events, the 
dimmer the stars they can affect. It is also to be noted 
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Table 5. Number of expected microlensing events for this 
pixel analysis and the EROS star monitoring analysis. The 
2nd row gives the number of events j\r^ APEROS expected 
given these efficiencies and a halo mass fraction /, that 
can be compared with the number of events iV™ os / / ex- 
pected on the same data (Renault 1996) shown in the 3rd 
row. The last row displays the ratio N£% APEROS /Nf£ os . 



M/Mq 0.01 0.05 0.1 0.5 1 

jyAGAPEROS/y 1Q2 Q77 Q 64 3 g 02g 

N™ os /f 0.14 0.055 0.045 0.018 0.014 

Ratio 6 11 13 19 19 



Table 6. Simulations for the EROS 91-94 data: the num- 
ber of events N^ APFjnos / f expected with a pixel analysis 
applying the same selection procedure as described above. 
The last row gives for comparison the number of expected 
events Nf^® s with the star monitoring. 



M 


0.1M© 


O.5M 


^AGAPEROS^ (pixel monitoring) 


5.0 


3.0 


N^t° S 1 f (star monitoring) 


0.42 


0.2 



6.1. The whole EROS-1 CCD data set 



that events affecting bright stars can be missed when the 
event duration is longer than the period of observation. 

5.2. Comparison with star monitoring 

Tab. |^ gives the number of microlensing events expected 
with our AGAPEROS pixel analysis as well as the num- 
ber of microlensing events expected with the EROS star 
monitoring analysis. For long-duration events due to M > 
O.5M Machos in the mass range where our analysis has 
been optimized, the number of microlensing events that 
our pixel analysis could detect is enhanced by a factor 
larger than 15. This gain is due to the fact that the EROS 
star monitoring analysis accounts for stars down to mag- 
nitude 19.5 but is far from complete down to the limiting 
magnitude, given the crowding and seeing conditions. In 
the LMC bar fields, stellar photometry has a bad detec- 
tion efficiency for dim stars, whereas our approach does 
not require to resolve the stars to detect their variation. 

Similar considerations can be made with respect to the 
MACHO star monitoring analysis (Alcock et al. 1996) ap- 
plied on a field 60 times larger than the 0.25 deg 2 field 
analysed here. The MACHO exposure was Emacho = 
9.7 x 10 6 star-yr. corresponding to an experiment du- 
ration of 409 days. In first approximation, the equiva- 
lent MACHO exposure for 0.25 deg 2 and 120 days would 
be E' MACHO — 4.7 x 10 4 star-yr. As the detection effi- 
ciencies present similar features as those discussed here, 
this exposure can be compared with the AGAPEROS one 
E agaperos = ^ s A £ APEROS x T ohs (yr) = 7.2 x 10 5 star- 
yr, which is about 15 times larger. 

6. Perspectives for the Magellanic Clouds 

We discuss here the perspectives than can be expected 
with the whole EROS-1 CCD data setJSect. |J), then 



with other data characteristics (Sect. 6.2), and give gen- 
eral considerations about further applications of this tech- 



The limitation of this analysis due to the short-time span 
over which the 91-92 data stretched could be overcome 
by analysing the whole EROS-1 CCD data set. We dis- 
cuss here the perspectives offered by such an analysis. We 
use the same simulations, but with the characteristics of 
the 6871 images available in blue and 7011 images in red, 
namely the time, absorption factor and background flux. 
Tab. ^ shows the number of events expected assuming that 
the same selection procedure^ can be used for the whole 
EROS-1 CCD data set. Hence, we can reasonably expect 
~ 4 x / events for 0.1 — 0.5Ajf Q Machos. This sensitivity is 
equivalent to the one typically achieved using a star moni- 
toring analysis performed on a much larger field (EROS-1 
plates: 25 deg 2 - around 10x/ expected events for 0.1 M© 
Machos). 

6.2. Other sources of data 

The pixel analysis described in this paper has been applied 
on data which were rather peculiar among the microlens- 
ing databases. An average of 10 exposures per night was 
available: this is the reason why we succeed to achieve 
a relative stability between 1 and 2% on super-pixels (see 
Paper I). For comparison, the typical level of photon noise 
obtained on the background of the EROS-2 images on the 
LMC Bar is about 0.3% (for super-pixels). Then the sta- 
bility that could be achieved on corresponding super-pixel 
light curves could be typically twice the photon noise, 
namely about 0.6%. Similar stability could be expected 
with MACHO images (D. Bennett, private communica- 
tion) . 

One must also keep in mind that this analysis has been 
applied to a field in the Bar of the LMC, which is a very 
crowded field. The gain of pixel monitoring would be sub- 
stantially smaller when used in less dense regions. 



nique (Sect. 6.3). 



3 In practice, the thresholds may have to be adjusted to ac- 
count for possible unexpected sources of noise, hence changing 
the sensitivity. 
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6.3. Pixel Method 

The principal problem with a pixel analysis towards the 
LMC is that, by its very nature, it has difficulty measuring 
the flux of the un-magnified star and the maximum mag- 
nification. One significant consequence of this problem is 
that it produces degeneracies that affect the determina- 
tion of the duration tg of the events. 

However, the same problem arises with star monitoring 
which is affected in a major way by blending: that is, a 
parameterisation similar to the one presented in Eq. [l] 
must also be considered (see Wozniak & Paczyhski 1997, 
Alcock et al. 1996) to account for the magnification of 
underlying stars. 

The MACHO and EROS groups (Alcock et al. 1997a, 
Pratt 1997, Renault et al. 1998, Palanque-Delabrouille et 
al. 1998) has corrected their events (detected by star mon- 
itoring) for blending effects with a statistical correction. 
Moreover, one could in principle overcome this difficulty 
with a high resolution image achieving a good signal to 
noise ratio down to V = 24 (see Han 1997). Ardeberg et 
al. (1997) have measured the flux of the stars in the LMC 
Bar on HST images down to magnitude 24 (Stromgren 
photometry) and claimed to be completed down to mag- 
nitude 22. Such a performance should solve the problem 
of measuring the flux of the un-magnified star for most 
of the events (see Fig. 14). However, the identification of 
the star that has varied is one possible practical prob- 
lem which must be studied further, but could probably 
be overcome. For the dimmest stars, unambiguous deter- 
mination of the stellar flux will most probably require an 
HST measurement of the star flux during the event. Such 
a requirement seems reasonable for an ambitious pixel ex- 
periment towards the LMC: HST measurements are al- 
ready being performed to correct for blending in events 
detected by star monitoring toward the bulge. 



7. Conclusions 

While the EROS-1 CCD data set has already excluded the 
small-mass Macho range, we have shown that it is possi- 
ble with a pixel analysis of the same data set to probe the 
mass range (0.01 — l.M©) where all the known events have 
been detected. We thus demonstrate for the first time the 
efficacy of the pixel method for LMC images. We have 
shown with the computation of detection efficiencies that 
the gain in detectable microlensing events with pixel mon- 
itoring is significant: pixel monitoring can detect 15 times 
more microlensing events than star monitoring for very 
crowded LMC bar fields. 

With simple selection criteria we have been able to 
show that no microlensing events are present in this data 
set. The criteria imposed mainly rely on the uniqueness 
of the variation, on the achromaticity and on the fact 
that the microlensing phenomenon - a geometrical effect - 
should affect all kind of stars independently of their type. 



This demonstrates that the noise which affects super-pixel 
light curves can be reduced to a level that is adequate to 
conduct the analysis. Our analysis, optimised for the de- 
tection of long duration events, has only detected variable 
stars, and we have shown our ability to reject such variable 
stars with the CMD. We have then demonstrated the ef- 
ficacy of this approach for the detection of variable stars, 
which remain the main background for microlensing. A 
comprehensive analysis of the variable stars rejected by 
this analysis is in progress (Melchior et al. 1998b). 

Such an analysis can detect microlensing magnifica- 
tion of unresolved stars up to V = 24 and is thus com- 
plementary to star monitoring. On-line pixel analysis and 
follow-up will probably be necessary to better discrimi- 
nate possible variable stars and to achieve a higher quality 
photometry. 

The short period of observation analysed here - 120 
days to be compared to 3 years for EROS-1 plates and 2 
years for MACHO - and the relatively small field - 100 
times smaller than the one analysed by the EROS group 
(plates) (Ansari et al. 1996) and 60 times smaller than the 
field analysed by the MACHO group (Alcock et al. 1997a) 
- explain the relatively small number of events expected. 
For the first time, we are able to provide detection effi- 
ciencies for a pixel analysis. These results allow us to esti- 
mate that a pixel analysis of the complete existing EROS-1 
CCD data base could detect about ~ 4. x / microlensing 
events in the mass range of interest (0.1 — 0.5M Q )... with 
a field of only 0.39 deg 2 ! 
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